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Abstract 

We derive a classical Schròdinger type equation from the phase space Li- 
ouville equation vàlid for an arbitrary nonlinear potencial V(x). The deriva- 
tion is based on a Wigner type Fourier transform of the classical phase space 
probability distribution, which depends on a small constant a with dimension 
of action, that will be characteristic of the microscopic phenomena. In order 
to obtain the Schròdinger type equation, two requirements are necessary: 1) 
It is assumed that the appropriately defined classical probability amplitude 
ty(x, t) can be expanded in a complete set of functions $ n (x) defined in the 
configuration space; 2) the classical phase space distribution W(x,p, t) obeys 
the Liouville equation and is a real function of the position, the momentum 
and the time. 

We show that for a equal to the Planck's constant h, the evolution 
equation for the classical probability amplitude ^/(x, t) is identical to the 
Schròdinger equation. The wave particle duality principle, however, is not 
introduced in our paper. 
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1 Introduction 



In 1932 E. Wigner published an important paper PjJ |2] where he intro- 
duced what is called today the Wigner's quasi probability function, or simply 
Wigner's phase space function. We shall denote it by Q(x,p,t) and it is 
connected with the Schròdinger wave equation solutions ip(x,t) by the ex- 
pression: 

Q(x,p,t) = — tp*(x - y,t)ip(x + y,t)e~^dy, (1) 
irn J 

where % is the Planck's constant. We recali that Wigner's intention was to 
obtain a quantum mechanical description of the classical phase space phe- 
nomena. 

It was stated later, by Wigner and O'Connel that the definition (0) is 
the only one that satisfies a set of conditions that are expected on physical 
grounds. 

Notice, however, that the Planck's constant h~ enters in the above expres- 
sion in a combination that was not completely clarified by Wigner and other 
authors [1-7]. We know that the Planck's constant h~ is related to the electro- 
magnètic fluctuation phenomena characteristic of Quantum Electrodynamics 
[H] and Classical Stochastic Electrodynamics [9-12]. However, for simplicity, 
these electromagnètic fluctuations will be not considered in our paper. 

We would like to mention that, if applied to an arbitrary solution ip(x, t) 
of the Schròdinger wave equation, the expression (JTJ) may lead to negative, 
or even complex, phase space probability densities, suggesting that many 
solutions of the Schròdinger wave equation cannot be associated with genuine 
phase space probability distributions, that is, the relation (0) is problemàtic. 
Simple examples are the excited states of the harmònic oscillator Here, 
following França and Marshall [H] , we shall consider that these excited states 
are simply useful mathematical functions. We shall clarify this point later. 

It is also well known [TJ|ü] that Q(x,p, t) defined in (JTJ) propagates in time 
according to the quantum Liouville equation 

dQ = _^dQ + dVdQ + (jkf&V&Q + {^d^Vd^Q + 
dt M dx dx dp 3! dx 3 dp 3 5! dx 5 dp 5 

provided that ip(x,t) satisfies the Schròdinger wave equation. Here F(x) = 
—dV/ dx is the arbitrary nonlinear force acting on the particle with mass M. 
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The equation (J2J) shows that the definition (JTJ), and the Schròdinger equa- 
tion, do not lead to the classical Liouville equation. We recali that, according 
to the correspondence principle, one should necessarily obtain the Liouville 
equation if the mass M is large enough. Therefore, the relationship between 
the equation ((TJ), the Schròdinger equation and the Liouville equation is also 
problemàtic. 

In our paper, starting with the classical Liouville equation, we shall see 
the mathematical and the physical conditions to connect it with an equa- 
tion which is formally identical to the Schròdinger equation for a classical 
probability amplitude t) defined in the configuration space. 

By means of the theorem established within the section 2, we show that 
the connection between the Liouville equation and the Schròdinger type equa- 
tion is vàlid for an arbitrary potential V(x). This connection was presented 
in reference [3], restricted to the case where the potential V(x) is quadratic 
in the variable x. 

Our presentation is organized as follows. We give, within section 2, 
the mathematical background necessary to connect the classical phase space 
probability distribution to the classical configuration space probability am- 
plitude t). For clarity sake, we shall denote ip(x,t) a solution of the 
Schròdinger wave equation and t) the classical probability amplitude. 
The direct relation between the Liouville equation for an arbitrary potential 
V(x), and the classical Schròdinger type equation is presented within section 
3. A brief discussion and our conclusions are presented in the final section. 

2 Mathematical background 

Our starting point is the Liouville equation for the probability distribution 
in phase space, denoted by W(x,p, t). The Liouville equation is such that 



The nonlinear force acting on the particle with mass M will be denoted 
by F(x) = —dV(x)/dx and V(x) is an arbitrary continuous function. The 
solutions of the equation Q are obtained in conjunction with the solutions 
of the Newton equation 






(4) 
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and the classical positive distribution function, W (x ,po), associated with 
the initial conditions Xq and po, namely 



W(x,p,t) = / dx / dpoWo(x ,po)6(x-x B (t))6(p-p e (t)). (5) 



Here x c (t) is the classical trajectory and p c (t) = mx c (t). Notice that 
x c (t), and p c (t), depend on x and po- Therefore, W(x,p, t) is always real 
and positive, and the classical probability amplitude *&(x, t) is such that 



The above definition is not enough to determine the complex amplitude 
ty(x, t). To achieve this goal we shall use the Liouville equation (J3J). 

The Liouville equation for W(x,p,t) has a close relation with a classical 
Schròdinger type equation for ty(x, t), as we shall see in the next section. In 
order to clearly explain this relationship we shall establish a useful theorem 
based on a convenient Fourier transform of W(x,p, t) similar to that used in 
equation (JTJ). 

Theorem : 

Consider the Fourier transform defined by 



where a is a small constant with dimension of action. The mathematical 
definition (J7|) has a physical meaning only when y pa because in this 
case, the equation (J2J) becomes identical to (jüj) and gives us the probabil- 
ity distribution in the configuration space. Mathematically, we can write 
T[W](x,y = 0,í) = |#0,í)| 2 . We shall show that T[W](x,y,t) can always 
be expressed in the form 




(6) 





T[W] {x, y, t) = V\x- y, t)*(x + y,t), 



(8) 



provided that W(x,p,t) is a real function of its variables. 



Proof: 
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Consider any set of functions <& n (x), that satisfy the completeness relation 
J2K(^n(y) = S(x-y), (9) 

n 

and are orthogonal to each other. An example of these functions is the set 



n{x) = ~JTb~' (10) 



where n = 0, ±1, ±2, ±3, ... and b is a positive constant (—b<x<b). Other 
sets of functions Q n (x) can be used without loss of generality [21 E]. 

Therefore, any complex probability amplitude ^(x,t) can be expressed 

as 

V{x,t) = 5> n (í)$„(x), (11) 

n 

where a n (t) are the coefficients of the expansion. Different sets of coefficients 
a n (t) will give different probability amplitudes ^(x, t). One can show, from 
the orthogonality properties of the functions $ n (x), that 

a n (0) = J dx$* n {x)y(x,0). (12) 

Let W mn (x,p) be phase-space functions defined by [21 E] 

W mn {x,p) = — í dy$* m (x - y)$ n (x + y)e~^ , (13) 

7TOí J 

where a is the same constant introduced in (Jlj). This is a convenient math- 
ematical defmition because the functions W mn (x,p) have useful properties 
as we shall see in what follows. Notice that the functions W mn (x,p) can be 
positive, negative or even complex functions [21 EI]- 

The functions W mn (x,p) constitute a complete and orthogonal set of func- 
tions in phase space. The completeness can be verified as follows: 

YlY, W mn(x,P)W m n(y,q) = 
m n 

Ot ) J J rm m 



[■na 



—5{x-y)5(q-p). (14) 
Tia 



In the last equality the fact that functions $ n (x) form a complete set was 
used, and we also used the Fourier integral representation for the Dirac's 
delta function. 

A verification of the orthogonality requirement can be obtained in an 
analogous way and leads us to the following relation E] : 

/ dx í dpW^ n (x,p)W rs (x,p) = —S mr 5 ns . (15) 
J J na 

Therefore, it is always possible to put the classical phase-space distribu- 
tion W(x,p, t) in the form 

W(x, P ,t) = J2J2 C mn(t)W mn (x, P ), (16) 

m n 

where the coeficients C mn (t) are only functions of the time t [ü]. 

By substituting the explicit form of W mn (x,p) in the last equation, we 
obtain the formal expression 

W(x,p,t) = — ídyY / T, C mn(t)K^-y)^n(x + y)e- 2 -^. (17) 

TTCï J m n 

Since W(x,p, t) is a real function we have 

W(x,p,t) = W*(x,p,t). (18) 
Using the relation (fT%J) we get the following equality: 

/ dyY.Y. C mn (t)$* m (x - y)$ n {x + y)e~^ = 

m n 

= í dyY.Y,C* mn m m (x-y)K(x + y)e 2 -^- (19) 

m n 

Considering the change of variables y = — £, and exchanging the dummy 
indices m and n, we can write (J19|) in the form 

/ d yJ212 c ™n(t)$* m ( x -y)®n(x + y)e-^ = 

m n 

= Jdt^^C^{t)& m {x-t)* n (x + t)e-—. (20) 
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From this equation we conclude that 

Cmnit) = C* nm {t). (21) 

Therefore, the elements C mn (t) can be written in the following form: 

C mn (t) = a* m {t)a n {t), (22) 

where the functions a n (t) are only functions of the time t jH] . In order to make 
the connection with the probability amplitude ty(x,t) we shall assume that 
the functions a n (t) are the same functions introduced above in the equation 

Using (J23), CU) and © we obtain 



dpW(x,p, t)e 



-•l'U 



— [dy> 
ira J 



5^a n (í)$„(x + y') 



Y,a n (t)<ï> n (x + y) 



dpe 



(23) 



So, according to the defmition Q and the expression for the classical 
probability amplitude, we get 



T[W](x,y,t) = / dpW(x,p,t)e— =V*(x - y,t)V(x + y,t), (24) 
thus completing the demonstration of the theorem. 



Considerations concerning the Liouville equa- 
tion and its connection with the Schròdinger 
equation 



Phase space distribution functions, as those defined in ([13|h provide a frame- 
work for a reformulation of non-relativistic quantum mechanics (QM) in 
terms of classical concepts [U 121 El- 
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Moreover, it is a widespread belief that the connection between the Liou- 
ville equation and the Schròdinger type equation is only possible for quadratic 
potentials in the x variable (V(x,t) = a(t)x 2 + b{t)x + c(t), where a(t), b(t) 
and c(t) are arbitrary functions of time [3]). We shall explain, in this section, 
that the potential energy V{x) can be an arbitrary function of x. In order to 
achieve this goal we shall use a procedure discussed before by L. S. Olavo jü] 
and K. Dechoum, H. M. França and C. P. Malta [7]. We shall see that the 
theorem established in the previous section is a proof of a working hypothesis 
(see ©) used in these and many other works U\- The first work quoted 
in the reference [7] treats the Stern-Gerlach phenomenon. 

As before, W(x,p, t) is the classical phase space probability density asso- 
ciated with some physical system (see the equations © and (jU)). We shall 
assume that W(x,p, t) is normalized so that 

J dx J dpW{x,p,t) = 1. (25) 

The classical configuration space probability density will be denoted by 
P(x, t) given by 

P(x,t) = J dpW{x,p,t) = |^(x,í)| 2 . (26) 

where ^(x, t) is the classical probability amplitude (see (jüj). 

We are interested in obtaining a dynamical differential equation for the 
classical probability amplitude ^(x, t) based on the fact that W(x,p, t) obeys 
the classical Liouville equation (jïïj). 

To obtain this equation, we shall use the Fourier transform 

/oo 2ivv 
W(x jPl t)e—dp, (27) 
-oo 

defined previously in the equation ((TJ. 

We recali that the Fourier transform T[W](x,y,t) is a complex function 
which has a physical meaning only in the limit \y\ — > 0, that is, 

T[W}(x,y = 0,t) = \q(x,t)\ 2 . (28) 

Therefore, we shall consider the definition (}2Tj) only for very small vàlues 
of y. Our goal is to obtain the differential equation for ty(x,t), from the 
differential equation for T[W](x,y,t), vàlid when \y\ is very small. 
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Our first step is to consider the equation 



^T[W](x,y,t) 



30 dW 2ipy 

— — e o, dp. 

00 at 



Using the Liouville equation (jïïj) we obtain 



dT[W] 
dt 

One can show that the integral 

I 

is such that 



p dW w ,dW 

— h Fix) 

M dx v ' dp 



t \x)— — e a dp 
dp 



2ipy 

e~^~dp. 



l = F[x)^-\ W(x,p,t)e^ ]l dp. 

Oi J-00 



Here we have assumed that 



lim W(x,p,t) = 0. 

\p\— *oo 



We also see that 



2ipy 

pe a 



a d 2ipy 



Substituting (jÏÏH) and (j33J) in (|UD we get 



ia^T[W](x,y,t) 



-iaf d 2 
2M dydx 



2yF{x) 



(29) 



(30) 



(31) 



(32) 



(33) 



T[W](x,y,t). (34) 



Notice that, in accordance with the theorem proved within section 2, 
T[W](x, y, t) = V*(x - y ,t)*(x + y, t). 

It is convenient to use new variables, namely s = x — y and r = x + y, so 
that the equation ()34|) can be written as 

d 



ia—[V*(a,t)V(r,t)] 



—ta ) 
2M 



D 2 



2 



dr 2 ds 2 



(r - s)F 



r + s 



V*(s,t)V(r,t). (35) 
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According to (|28|). we want an equation for *&(x, t) = [\I/(r, t)} y __> , and 
the equivalent equation for t) = [^/*(s,t)] y ^ . Consequently we shall 

consider that the points r and s are arbitrarily close so that 



- (r - s)F 



£ F(OdÇ = V(r)-V(s), 



(36) 



in accordance to the mean value theorem. Substituting (jïïü|) into íjïïo} we get 
the equations 



and 



d 

ia— V*(s,t) 



ia— tf(r, t) 



2M 



J— \-ia^-\ + V(r) 
2M \ dr y 1 



V(r,t). 



(37) 



(38) 



These two equations are, in fact, the same differential equation, namely 

2 



ia—W(x,t) 



\ 1 ( 




2M \ 





V(x) 



(39) 



which is formally identical to the Schròdinger wave equation. Notice, how- 
ever, that fy(x,t) is a classical probability amplitude, which is conceptually 
different from the Schròdinger wave function ip(x,t). The classical probabil- 
ity amplitudes *f>(x,t) are just mathematical objects [Tü] . 

The wave-particle duality principle was not used and the small charac- 
teristic constant a is, up to this point, arbitrary. This constant can be 
determined from the observation of several phenomena within microscopic 
domain [10-12]. Therefore, we can conclude that a — H, where h is the 
universal Planck's constant. 

It is relevant to observe that if we write 

*(x,t) = $ n (x)e-^, (40) 
where e n are constants with dimension of energy, the equation (|39|) leads to 



2M 



-ia 



dx 



V{x) 



$n(x) = e n $Jx). 



(41) 
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This is a familiar eigenfunction equation . Its solutions $> n (x) constitute 
a complete set of orthogonal functions that can be obtained for different 
potentials V(x). This set of functions <& n (%) can also be used in our previous 
equations (fTTj) and (fTïïj) (see reference jÜ]) showing the consistency between 
the sections 2 and 3. 

We would like to call the reader attention to the fact that the above 
mathematical treatment can be extended to the three dimensional case [7] 
without additional difficulty. 

4 Discussion 

The connection between the classical Liouville equation for a probability 
distribution in phase space and the Schrodinger equation was already es- 
tablished for quadratic potentials 0] . The extension of the formalism to 
a genèric potential V(x) presented here put all the previous related works 
[7j on a sòlid mathematical ground. We get the Schrodinger type equation 
and also Born's statistic interpretation of \^f(x, t)\ 2 , by considering the 
Liouville equation, and the equations (J7J), (|34j) and (|35|) only in the limit 
y — > 0. In this context we were able to use the mean value theorem (see 
equation (jHüJ)) to separate (|33|) in two equivalent equations, formally iden- 
tical to the Schrodinger equation. Since y has to be small, we understand 
why the inverse transform (equation utilized by Wigner, does not guar- 
antee that the phase space function Q(x,p,t) is a true (positive) probability 
distribution. 

The association of the constant a with the Planck's constant % is natural 
and inevitable to any reader familiarized with the Schrodinger equation used 
in the microscopic world. As we said above, the explicit numerical value of 
a (or h) can be determined by testing the validity of the Schrodinger type 
equation in various phenomena of the microscopic domain [T^j. The connec- 
tion of the Schrodinger momentum operator —ihd/dx with the zero-point 
vacuum electromagnètic radiation is also interesting. This fact is explored 
by P. W. Milonni in the references [TJJ ^3 • 

According to A. J. Legget [Tïïj, "despite the spectacular success of quan- 
tum mechanics over the last 80 years in explaining phenomena observed at 
atòmic and subatòmic level, the conceptual status of the theory is still a 
tòpic of lively controversy" . He seems to believe that quantum mechanics is 
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nothing more than a very good mathematical "recipe" |17j . 

We agree with the above statements. According to the connection be- 
tween the Liouville equation Q and the Schrodinger type equation (JHHJ), 
presented here, one can conclude that the classical probability amplitude 
ty(x, t) is not associated with a genuine de Broglie wave. As a matter of fact, 
the existence of de Broglie waves was questioned recently by Sulcs, Gilbert 
and Osborne ^H] i n their analysis of the famous experiments ^ü] on the 
interference of massive partides as the 6*60 molècules (fullerenes). 
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